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A method is proposed for using sequentially-fired, shaped charges to drill 
deep, small-diameter holes in rock more efficiently. Using this method, 
penetration in igneous rock can be increased from 10 times the diameter 
of a single, copper-lined charge to 18 or more times this diameter when 
two or more charges are fired sequentially. Design calculations for a shaped- 
charge drilling system can be performed using Bernoulli's equation; however, 
this simple mathematical model has several shortcomings. Shaped-charge 
penetration, therefore, should be computed using a scaling technique and 
available experimental data. This allows the penetration efficiency of a 
sequence of shaped charges to be determined. A shaped-charge drilling device 
is presented which uses two or more staged charges. Hydrodynamic evalua- 
tion of this device, especially the mechanics of slug deflection, indicates that 
this design overcomes some of the weaknesses of earlier designs. A major 
improvement is the use of explosive cylindrical lenses, initiated by mild 
detonating delay fuses, to ensure proper symmetrical firing of the second- 


and third-stage explosives. 


INTRODUCTION 


Shaped-charge drilling, the most rapid drilling process 
known, is a spinoff from World War II technology de- 
velopment. A shaped charge is an explosive with a 
metal-lined, hemispherical-, wedge-, or conical-shaped 
cavity. Detonation collapses this cavity, ejecting a 
supersonic stream of metal called the jet [1]. Drilling 
speeds of 2-6km/sec and hole depths of 4-15 times 
charge diameter can be achieved using shaped charges. 
Despite this tremendous potential for mining applica- 
tions, shaped charges are used primarily in cutting and 
perforating metal. Their potential for tunneling and 
blasting work should be considered seriously, since 
such explosives provide a very concentrated energy 
source. The shaped charge lends itself to remote mining 
operations where noise is not an important factor. 
High-velocity penetration of solid materials, as 
opposed to slow mechanical drilling or ballistic 
penetration, involves the melting and vaporization of 
the target. As a result, penetration depth is little depen- 
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dent on the strength of the material being penetrated. 
Simple scaling laws, therefore, can be used quite effecti- 
vely in predicting relative jet effectiveness in penetrating 


‘different materials. 


THEORY 


Simple model 


The most conceptually useful description of target 
penetration by a shaped-charge jet (Fig. 1) and, for 
most purposes, a reasonably good one is based on Ber- 
noulli’s equation 


é 
<P + 4U? = constant ~ ee (1) 


where p is pressure, p is density, and U is velocity. 
The last term of this equation is based on the assump- 
tion that the target material compresses only slightly 
or not at all during the process. Then, following Birk- 
hoff et al. [1] the equation can be rewritten as 


3p({V — UP = 3pU", (2) 


where p is the density and U the velocity of the target 
material at the jet-target interface and V is initial vel- 
ocity of the jet with a density p, (see Fig. 1). Since 
penetration time (At) is the time it takes for the rear 
end of the jet to reach the moving interface, 


l 


At = V—U (3) 

and, using equation (3), we represent penetration by 
U ! p; 1/2 4 
P= UAt=I>—G = (2) : (4) 
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Ejected fragments and vapor 


Fig. |. Diagram of idealized jet penetration (V = jet velocity; U = 
penetration rate). 


Penetration depth, therefore, depends only on jet 
length. 


Approximations 


The penetration model just described is not com- 
pletely accurate for several reasons. One, the jets pro- 
duced from shaped charges do not have a constant vel- 
ocity along thei: length [2]. This means jet Jength 
changes with time or travel distance, so that the jet 
eventually breaks up into a number of elongated par- 
ticles, each acting as an individual penetrator. These 
particles, in turn, begin to tumble, thereby reducing 
their penetration capacity. Breakup and tumbling are 
the results of perturbing oscillations produced in jet 
formation and turbulent interactions with surrounding 


gasses, respectively. These factors are not taken into 


account in the simple penetraton model. 

The shape of the jet particles also affect jet length 
and, hence, penetration depth. Particle shape depends 
on the material used for the liner. X-ray shadowgraphs 
[3] indicate that an iron liner breaks up like a brittle 
material, forming jagged particles; while copper 
appears more ductile, breaking up into elongated par- 
ticles. Aluminum jets, on the other hand, appear like 
a spray of particles down to submillimeter size probably 
with smaller, unresolved particles in between. This sug- 
gests partial melting of the aluminum. 

Finally, the mode! ignores the transfer of momentum 
from material ejected from the target hole to the im- 
pinging jet. It also ignores viscosity effects at the target- 
jet interface. 


Jet length 


The fact that jet velocity is not constant and jet 
length varies with time and distance increases the value 
of | in equation (4) with time. The resulting increase 
in penetration capacity with increased standoff (dis- 
tance between charge and target) is well known [1]. 
Penetration—vs-standoff curves have the following 
general shape: an initial increase of penetration (P) as 
charge and target are separated, with maximum 
penetration occurring at separations between 7-10 
cone-base diameters, followed by a decrease in 
penetration with greater standoff distances. 

The performance of a precision-machined copper 
cone has been examined experimentally using several 
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Fig. 2. Shaped charge performance (copper cones and steel target): 
penetration vs charge to center of hole in cone-base diameters (- = 
Octol explosive; A = HMX-based extrudable explosive). 


different explosives. As shown in Fig. 2, Octol and a 
(60:40=) RDX: TNT mixture both yield jets with 
about the same steel penetration capacities [4]. Impor- 
tant to note is the decrease in efficiency as charge-to- 
target distance is increased. A similar effect was found 
using a more energetic charge of extrudable HMX-base 
explosive developed at Lawrence Livermore Labora- 
tory [5]. However, in this latter case, performance im- 
proved at nominal standoff distances and dropped off 
slightly faster at greater distances. The flags in Fig. 2 
indicate the range of penetrations obtained from three 
or more charges. This variability shows the importance 
of quality control. Ideally, it should be possible to pre- 
dict penetration into any other target by scaling these 
data with equation (4) but minor deviations may be 
expected in low-density targets. 


Parameterization 


A variety of parametric equations to account for jet 
lengthening and breakup [1,2], change of cross sec- 
tional areas [1], and target strength [6,7] have been 
formulated and compared with experimental data. 
Detailed X-ray analysis of jet penetration [8,9] has 
shown that jet length increases to the time of breakup 
and that subsequent particle length (hence, effective jet 
length) remains constant. The latter holds true up to 
the time that the jet particles start to tumble, reducing 
the effective value of |. A strength factor was introduced 
by an empirical determination of the jet particle vel- 
ocity below which no penetration occurs. It must be 
noted that this limiting velocity is not only a function 
of target strength, but also of target density and stand- 
off (as may be inferred from the above discussion 
and equation (4)). Penetration data indicate, however, 
that target strength is not an important variable in des- 
cribing variations between materials. 


Alternate approach 


A detailed analysis of the penetration [10] process 
has been made with limited success [9]. Actual 
penetration is considerably smaller than that predicted 
by this simple theory, which accounts for jet length 
variations and breakup, but not the other limitations 
of equation (4). A somewhat different approach to com- 
paring penetration through materials with different 
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densities is to rewrite equation (4) in the form 


dP = (pj/ps1)* OI (d), (9) 


recognizing the fact that efficiency is a function of dis- 
tance and that we wish to use a value of (=/*") consis- 
tent with observed penetration yet smaller than the 
actual value of |. It is convenient to make comparisons 
to the peak of the P/D-vs-d/D curve in Fig. 2; a numeri- 
cal method is employed to avoid fitting bell-shaped 
curves. 

Considering the penetration of n sections of jet, each 
of length di‘ (d), we get 

Pse 2, OP (di) = 0 (0)/ Ps)? . SI (d)). (6) 

Now, realizing that the value of d/*' is a maximum at 
some value of d=d,, we may for any real process 
calculate a maximum penetration (P™*) by normalizing 
each term by 


dI"(d1) _ P.Ady), 
dd) P,{d) 


a(d) = 


which can be calculated from the steel penetration curve 
(Fig. 2) with reasonable accuracy. Therefore, 
Pom = Y (pila) 6d) afd) (7) 
represents this maximum penetration in steel, which may 
be compared to the penetration through any material 
(m) by multiplying the terms by (p,,/p,,)'/2. We can now 


write the penetration-scaling equation between steel and 
material m as 


Pa = Y (Pi/Pm)'? -(Pm/Pa)”” - 5°%(di) . afd) 


ll 
rer 


Be ¢ OP,,{d;) 7 a(d;), (8) 
where, in the last sum, B,, = (Om/Py)'’? and SP,, is the 


incremental penetration of the material in question. 
é 


EXPERIMENTAL CHECK 


To test equation (8) in an extreme case, a precision 
charge like that shown in Fig. 2 was fired through a 
0-61-m column of water into a stack of steel plates. 
The cone-base diameter was 0-053 m. Penetration and 
jet configuration were monitored by X-ray shadow- 
graphs. In addition, penetration was monitored with 
a set of charged foil switches that recorded the position 
of the jet tip at various times. The experimental setup, 
including significant dimensions, is shown schemati- 
cally in Fig. 3. A simple device called the virtual origin 
approximation [9] allows us to calculate the position 
of the jet particles quickly as a function of time. This 
surprisingly accurate calculation assumes that all par- 
ticles start at the same time from a common origin. 
The progress of the jet particles through the water and 
steel layers is shown by the dashed lines in Fig. 4. The 
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Fig. 3. Schematic of penetration tests. 


position of the jet tip, obtained directly from foil switch 
closing times, also is shown. This position can be in- 
ferred, in part, from the X-ray shadowgraphs that are 
superimposed on the graph. These shadowgraphs corre- 
late well with the penetration locus given by the foil 
switches and show the shape of the particles. 

Now, equation (8) can be used to calculate maximum 
penetration into the steel from the above experiment 
for comparison with Fig. 2. Using B,, = 0-357 and tak- 
ing a from Fig. 2, we get p{°* = 7.00, in good agree- 
ment with Fig. 2. A word of caution must be added, 
however. In extrapolating penetration in steel to that 
in water or rocks, allowance must be made for the dif- 
ference in the fraction of the jet effective in penetration. 
In the above example, this cutoff was at a jet velocity 
of 4km/sec (see Fig. 4). But this holds for the 
penetration through steel. For soft, low-melting point 
materials the cutoff velocity will be lower and effective 
penetration higher, an important factor especially in 
soft or hydrated and water-bearing rocks. 

To estimate the effective penetration depth of a 
sequence of charges in material other than steel (e.g. 
rock of density 3 g/cm), we modified equation (8) and 
the graph in Fig. 2 by rewriting the left side of the 
equation as 


P 
bPmax such that 6Pmax —_ st 
i=1 ra 


Me. 


(9) 
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Fig. 4. Jet characteristics in water and steel; dashed line shows position of jet particles vs time in the virtual origin 
approximation; plus signs indicate foil switch signals produced by jet tip passage; dots represent jet particles. 


Equating terms, dividing through by x and f, and 
summing we get 


n n dpmax 
P,, = Y OP = (1. Bm) |. Ye, 10 
p> . B ) p> a(d;) a 
where 
| p-l Pmax Pm 
ai ai i 2a Tak In. ak 


A reasonably good estimate may be made with n = | 
or 2. The first of a sequence of 50mm-dia charges 
would penetrate 0:5 m, the second stage an additional 
0-4 m, and a third charge an additional 0-2 m; making 
total penetration for the two-stage device 09m or 18 
diameters and for the three-stage device I-1m or 22 
diameters. The limitation of multistage drilling devices 
clearly is the reduction of penetration as a function 
of distance. 


STAGING CONCEPT 


Hole size and drilling depth can be controlled to 
some extent by standoff and cone size. Increased stand- 
off decreases hole diameter but increases hole depth 
(up to a maximum distance of 10 cone diameters). In- 
creasing cone size increases both diameter and depth 
of the hole. For many applications, however, increased 
penetration is more important than increased diameter. 
Furthermore, where space or total explosive weight are 
limitations, staged charges (i.e. firing in sequence a row 
of carefully aligned charges) should be considered [4]. 


A number of staged, shaped-charge drilling tech- 
niques have been proposed. See, for example, U.S. 
patent: 3,215,074 (2 Nov. 1965); 3,375,108 (26 March 
1968); and 3,416,449 (17 Dec. 1968). The multistage 
drilling device shown in Fig. 4 uses a normal conical- 
shaped charge (A) with a 42° apex angle and as many 
charges with hollow truncated cones (B) as the appli- 
cation and ultimate design require. The firing sequence 
is as follows: charge A is detonated by a single detona- 
tor centered on the exposed end of the explosive charge 
(not shown). This detonation collapses the metal cone, 
which forms a jet that passes through the central tube 
of the slug catcher between charges A and B. The 
detonation wave from charge A ignites an explosive 
delay train, which lights an explosive charge in the slug 
catcher, trapping the slow part of the jet after the main 
portion of the jet has passed. Meanwhile, the explosive 
fuse train burns on to light another explosive delay 
coil that, in turn, ignites a cylindrical lens which in- 
itiates a symmetrical detonation wave in the second 
shaped charge B. The process, in principle, may be 
repeated for additional charges. 

The device can be sealed and evacuated to eliminate 
the perturbing effects of air shocks; however, there are 
some problems yet to be solved. These include: 


1. Removal of the slow part of the jet (slug catching); 

2. Mechanical isolation of charges to prevent un- 
wanted shock deformations and preinitiation of neigh- 
boring charges; 

3. Proper initiation of the sequenced charges. 
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Shaped charge Slug catcher 
A) 


Shaped charge 
(B) 


Fig. 5. Schematic of multistage, shaped-charge drilling device (e = 
explosive; Cu = copper cone; p = porous metal shock attenuator; 
d = mild-detonation-fuse delay coil; | = cylindrical explosive lens). 


SLUG CATCHER PROBLEM 


The slug catcher design shown in Fig. 5 is schematic 
only. Neither the central tube, the explosive, nor the 
shape of other downstream parts need to be cylindri- 
cally symmetric. Furthermore, few materials of practi- 
cal interest exist that can be penetrated by copper jet 
particles moving at 2km/sec. In fact, at this low vel- 
ocity the trailing metal (slug) tends to fill the hole made 
by the fast portion of the jet. The slug catcher, there: 
fore, must be designed to trap a high-density metal 
stream flowing at up to 4km/sec, recognizing the fact 
that much of the metal cone (~ 60-70%) is accelerated 
to velocities between 1 and 4km/sec. 

If the slow portion of the jet is to be trapped by 
symmetrically compressing the tube, the energy density 
and momentum of that part of the jet must be consi- 
dered very carefully. Copper jet speeds of 3-5 km/sec 
or greater are required to penetrate steel. Slug trapping 
by symmetric tube closure around a jet depends on 
the viscous transfer of momentum to the tube from the 
jet stream either close to its melting point or in the 
liquid phase. This viscous deceleration process requires 
a relatively large amount of explosive to maintain tube 
pressure during the deceleration period. One way to 
increase the deceleration rate is to shape the charge 
to give material surrounding the jet a velocity com- 
ponent in the opposite direction of the jet. This method, 
however, requires a comparatively large amount of 
explosive. In addition, the design complexities involved 
could double system costs. 

A better approach is to introduce a small amount 
of momentum perpendicular to the jet to deflect it 
against a sufficiently heavy mass. A simple way to cause 
such deflection is to detonate asymmetrically a cylinder 
of explosive around the steel tube through which the 
jet passes towards the second stage (see Fig. 5). This 
design technique was tested by calculation and exper- 
iment. The calculations were performed using a 
lagrange hydrodynamic code [11], which followed the 
behavior of all significant parts of the system as the 
explosive detonated, compressed, and accelerated sur- 
rounding parts. This calculated behavior at the start 
of detonation and when the hole is almost closed is 
shown in Fig. 6. The deflection of the tube from its 
axis is small for the design dimensions chosen. This 
small deflection was also confirmed by experiment. In 
fact, recovered samples suggest that later in time tube 
deflection decreases slightly. 

A more effective design is depicted in cross-section 
in Fig. 7. This is a calculation of the movement of 
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Fig. 6. Calculated implosion of a solid iron pipe using asymmetric 
explosive initiation. 


a slotted cylindrical system of explosive and steel. The 
particular configuration chosen for the calculation 
results in the central jet being enveloped by the tube. 
The choice of slot angle or shape is clearly in question 
here and depends to some extent on the overall design. 
In any case, the central jet material reached a radial 
velocity of 1 km/sec, 5:5 usec after detonating the 
catcher explosive. This is more than enough velocity 


. Calculated jet deflection by a slotted catcher system. 


High explosive 


Jet path 


Fig. 8. Slotted slug catcher configuration. 


for jet deflection. This shape, therefore, appears to be 
a good guide for the design. 

An important consideration is the effect that severing 
the rear part of the jet could have on the remaining 
portion. Some of the jet material will be partially de- 
flected. The design challenge is to minimize this effect. 

A second consideration is the strength of the catcher 
material after deformation by the explosive. This is im- 
portant because it must withstand the impact of mas- 
sive, slow-moving slug fragments that could reopen the 
jet flight path and exposing the second stage to slow 
debris. Final slug catcher design must account for both 
of these possibilities. A reasonably simple, but ade- 
quate, design of a slotted catcher system is shown in 
Fig. 8. 


SECOND-STAGE INITIATION 


If protected from debris, the second stage can be 
detonated symmetrically and in proper alignment with 
the first stage. This second-stage initiation can be 
accomplished in two ways. If limited space is not an 
important factor and a high-voltage source is readily 
available, the second stage can be initiated using a ring 
of detonators. This method, however, is bulky and cum- 
bersome to use under field conditions. The second 
method is to use mild detonating fuse between the two 
stages and an explosive lens to ignite the second stage 
symmetrically at {0 or more points on a ring around 
the charge. In addition to being more compact, the 
cost of handling and deploying this system should be 
lower. Moreover, since the lens consists of a system 
of grooves cast in a plastic cylinder and filled with 
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extrudable or plastic explosive, it readily lends itself 
to the economies of mass production. 


SUMMARY 


This study shows that the. drilling performance of 
shaped charges can be improved significantly by stag- 
ing the charges and firing them in sequence. Based on 
hydrodynamic calculations and actual experimental 
tests, a multistage, shaped-charge drilling device is pro- 
posed that satisfies the design restrictions imposed by 
metal jet behavior. These design calculations also 
demonstrate the value of hydrodynamic computer 
codes in designing drilling systems of this type. Finally, 
an analytical method is developed and evaluated, which 
allows shaped-charge drilling performance to be scaled 
to a wide range of target materials. 
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